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The amount of screening of a proton in a metal, migrating under the influence of an 
applied electric field, is calculated using different theoretical formulations. First the lowest 
order screening expression derived by Sham» is evaluated. In addition 'exact' expressions 
are evaluated which were derived according to different approaches. For a proton in a metal 
modeled as a jellium the screening appears to be 15 ± 10%, which is neither negligible not 
reconcilable with the controversial full-screening point of view of Bosvieux and Friedel^. In 
reconsidering the theory of electromigration, a new simplified linear-response expression for 
the driving force is shown to lead to essentially the same result as found by Sorbcllc 3 , who 
has used a rather complicated technique. The expressions allow for a reduction such that 
only the scattering phase shifts of the migrating impurity are required. 

Finally it is shown that the starting formula for the driving force of Bosvieux and Friedel 
leads exactly to the zero-temperature limit of well-established linear response descriptions, 
by which the sting of the controversy has been removed. 

I. INTRODUCTION 

The amount of screening of a hydrogen atom in a current carrying metal has been the subject of a long-standing 
controversy. In brief, considering the driving force F on such an atom as being composed of two contributions, a 
direct force Fd due to the charge of the proton and a wind force F w due to transfer of momentum of the current 
carrying electrons to the proton, so 

F = F d + F w = (Z d + Z w ) eE = Z*eE, (1) 

Bosvieux and Friedel 2 found a complete cancellation of Fd, implying full sceening of the proton charge and only a 
wind force being operative, while most other researchers in the field were in favor of at most a very limited screening^ 
According to the convention in electromigration theory, the forces in Eq. are written as being proportional to 
corresponding valences and the applied electric field E. The effective valence Z* is the measurable quantity. The 
wind force has been calculated reliably for many systems with ab initio methods for the electronic structure. This has 
been done not only for migration of interstitials such as hydrogen, but also for substitutional impurities, including 
self-electromigrationiii£i£iL& 

In view of the long history of the controversy regarding the direct force we mention just a few key papers. The linear 
response expression for the driving force derived by Kumar and SorbelloS was considered as a sound starting point 
for the resolution of the controversy. From an evaluation to lowest order in the impurity potential Sham concluded 
to a negligible screening. 1 Using an evaluation up to all orders in the potential, in 1985 Sorbello found a screening of 
at most 25%£ This result, combined with a measurement of negligible screening in V(H) and Ta(H) and a screening 
of the order of 50% in Nb(H), led to a consensus^ Nevertheless, Friedel kept defending that only a wind force was 
operative, the more so as Turban et al. had given another support for that point of view— The confusing feature of 
the latter work is that their starting formula is a well-established form of a linear response expression. However, these 
authors do not evaluate that expression, but they use a proportionality argument regarding the expression for another 
physical quantity. Another support for the full-screening point of view was given by the present author—. However, 
that result was considered to be valid in the low temperature limit onljii&i^, and this limit is a rather academic one in 
view of the relatively high temperatures at which electromigration experiments are carried out. More recently Ishida 
predicted a screening ranging from to 100%, but his results were depending sensitively on the electron density of 
the host metal— 

We will present a thorough study of the amount of screening. After summarizing the main ingredients of the 
linear response description of the driving force in SjH] this is done by first evaluating Sham's screening expression 
numerically for a number of model potentials representing the impurity. The results, given in £11111 do not support 
Sham's conclusion of a negligible screening, but they are in line with Sorbello's result of a screening of 10-30%. 

Secondly, in illVI we will present a very simple evaluation of the linear response description. This evaluation is 
supplementary to the evaluation given by Rimbey and Sorbello^A^ and furthermore much more straightforward. The 
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two descriptions are compared in SjV] In flVIl it appears to be possible to reduce the final expression for the direct 
force valency Z d to a form containing just the scattering phase shifts of the migrating impurity potential. Numerical 
results will be presented in Will and compared with Sorbello's results^ 

Finally, in Willi we will show that the starting expression of Bosvieux and Friedel for the driving force is precisely 
the zero temperature limit of well-established linear response expressions. This is found by describing the switch-on 
of the electric field properly and by giving credit to the hermitian property of the Hamiltonian of the unperturbed 
system. 



II. LINEAR- RESPONSE DESCRIPTION 



The linear response expression for the force on an impurity with chemical valency Zi at a position Ri due to an 
applied electric field is given by 

F = Z t eE - ieE v dt e - at Tr{p(H) F op (t), £ } = F^' + F^ otal . (2) 

The first term clearly is the direct force on the bare ion. The cartesian label v runs from 1 to 3, the infinitisimally 
positive number a represents the adiabatical switch-on of the electric field represented by the potential 

8V(t) = eEe Q * • ( ^ rj - Z t ^ R Q ) = SVe at , (3) 



with j running over the electrons and a over the ions, and the operator p(H) is the grandcanonical density depending 
on the system Hamiltonian H . The force operator contains the electron-impurity potential 

V Bi = Y d v{T j --R a ) = Y t v*, (4) 

which is part of the system Hamiltonian, and is given by 

F p = -V Rl V ei =-^V Rl t)( rj - RO =J2 f j- (5) 

3 j 

Its time dependence refers to the Heisenberg representation 

F op (t) = e iHt F op e- iHt . (6) 

It appears that the second term in Eq. which is of course supposed to lead to the wind force, also contains 
some screening contribution to the direct force. The controversy has not to do with the fact that there is a screening 
contribution in F^ tal , but it is as to the magnitude of that screening contribution that people don't agree. The 
expression published by Kumar and SorbelloS, 



jrtotal = _i Eu I dte - 



at Tr{p{H) 



(7) 



follows simply and straightforwardly from a partial integration of Eq. with respect to the time. The current vector 
is defined as 



J=ie[5>,2?] = -e£^ = £^ (8) 



The driving force (J2J can be decomposed as follows: 

F = Z t eE + F^ otal = Z,eE + F; cr + F|f = {Z t + Z scr + Z w )eE = (Z d + Z w )eE, (9) 

containing the result of Bosvieux and Friedel for the wind force F^ F and a screening contribution 2 . In all treatments 
available F^ F can be written in its general form 



fS f 



- J 5n(r) r; 7B. 1 v 1 d 3 r, with u 1 =u(r-R 1 ). (10) 
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The precise explicit form depends on the level of approximation used to represent Sn(r), which is the local deviation 
of the electron density from its unperturbed host value due to the applied field and the presence of the impurity. 
From now on we will concentrate on Z c \. 

All previous relevant descriptions have been given for the electron-impurity system, for which the Hamiltonian H 
can be written as a sum of single particle Hamiltonians h, so 

H = J2 h3 with h = h +v = h + J2 ya - ( n ) 

3 a 

This allows for a reduction of the many body expression in Eq. to the following single particle expression, 

Fw' al - -ieE v J°°dt e- at tr{ [r v , n(h)]f\t)}, (12) 
where n{h) is the Fermi-Dirac distribution function in operator form 

It has been shown explicitly that if in the right hand side of Eq. 11211 the statistical operator is replaced by this 
operator for the free particle system, so n(h) — > n(ho), the Bosvieux-Friedel wind force expression F^ F arises 17 . That 
means that the screening part is given by 

¥ s ™ = -ieE v J^ dt e- at tr{[r v ,n(h) - n(/i )]f = Z scr eE. (14) 

The screening valency Z SCI is defined as 

Z SCI = -1^"^ e- at tr{[r,n(h) - n(h )] -f\t)}, (15) 

in which the factor of | comes from the fact that all three terms in the inner product of the vectors r and f 1 contribute 
equally. In all further evaluations the metallic host is modeled by a jellium, which is the only model used so far in 
the literature for the study of the direct force problem. This means that the electrons are perturbed by the random 
distribution of impurities only. Following Shami we now first consider the result to lowest (second) order in the 
impurity potential v. 



III. EVALUATION OF SHAM'S EXPRESSION 

The evaluation of Eq. Ijl5|l to lowest order in v requires the expansion of the statistical operator n(h) — n(ho) in v, 

n(h) = n(h ) - n(h) ds e sh v e - sh °(l - n{h )), (16) 
Jo 

while in the time dependence of f 1 one can replace h by ho. One obtains 

zscr = 4 E _ k l^ j* / dnu _ 

3m ^ (e k - e k >Y + a 2 \de k e k - e k , J 

where k is a free electron wave vector and e k is the corresponding energy. The matrix element < k|[r ly ,n(/i) — 
n(/io)]|k' > is most easily evaluated if one realizes, that it is equal to «(gf- + gfr) < k\n(h) — n(/i )|k' >. Following 

Sham and Sorbello^ the potential v refers to the migrating impurity only. Sham stored part of the presence of the 
impurities through the relacement a — * r _1 , t being the transport relaxation time due to the impurities, which can 
be justified by an average over the distribution of the impurities in the time dependence of the force operator. Both 
Sham and Sorbello were able to make their complete derivations after taking the T — > limit only. It has been shown 
that Eq. (|17fl reduces to Sham's expression after taking that limit ^ 

A numerical evaluation of Z scr becomes possible if one employs the spherical wave expansion for a plane wave, 
converts the summations over the wave vectors to integrals and carries out the angular integrals over the directions 
of the wave vectors. After using the relation between k 2 and the energy e k one ends up at 

Z- = - / de k de k/J ^ ^±_STMk,k'), (18) 
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in which the function ft(k, k') is defined as 

fi(k, k') = e k y/tvvi(k', k) 



(21 + l)kv e (k', k) - 2(1 + l)k'v e+1 (k', k) 



containing the information about the ion potential through 



Vi(k', k) 



r 2 dr jt(k'r)v(r)ji(kr) 



(19) 



(20) 



The integrand has to be treated with care when e k > = e k , because then the denominator attains the value a 2 which 
would imply 'singular' behaviour. However, precisely then the numerator becomes zero, because lim^ ,_» £fc (n k — 
n k >)/(e k ~ £ k') — * The crucial part of the integrand lies in the square around the point (e k ,e k ') — (ef> £f)- hi 

studying the Z scr integral it appears that in that square one has to keep the Fermi-Dirac distribution function in its 
finite temperature form. We could not obtain a reliable stable numerical result by using Sham's T — > expression. 
The result of a numerical evaluation for different ion potentials is shown in FIG. ^ We used a screened Coulomb 




FIG. 1: The amount of screening represented by Z scr according to Eq. I|18|l . for the screened Coulomb potential and for two 
square well potentials. 



potential 

v ( r ) = - Z ' e2e ^ while v k , k = -L fd 3 r e l(k - W > r v(r), (21) 
with Zi = 1 representing a proton in a jellium and the inverse screening length A given by the Thomas-Fermi 

1 9 

expression 

A 2 =W^(e F ) with M(e) = 7 ^^- -> v kk = ~ Q -g. , . (22) 

In addition square well potentials were employed in the same spirit as Sorbello did£ The width ro of the square 
well potential was chosen to be equal to the screening length 1/A and twice as large. The corresponding well depth 
vo was determined by the condition vo/v c = 0.999, where v c = 7r 2 /8mr 2 is the critical value of the well depth for 
which a bound state forms. For further details, see RefA The value of A is determined by the Fermi energy. While 
Sorbello chose five values for the Fermi energy, typical for metals ranging from sodium to aluminum, we have done 
the calculation for a whole range of Fermi energies. The results are plotted as a function of the Fermi wave number 
kp. The fcp values of sodium and aluminum are indicated. 
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FIG. 2: v\ h plotted as a function of /cf, for the two square well potentials and the screened Coulomb potential. 

Because A increases monotonically with the Fermi energy, the range of the corresponding screened Coulomb potential 
decreases with increasing fcp, which reduction in strength is seen clearly in the solid curve. In FIG.[21the v\ k curves 
are plotted for the three potentials. A clear decrease is seen for the Coulomb potential, and a rather flat behaviour 
for the square well potentials, while the one with 2/ A is markedly stronger than the one with the smaller width. 
Apparently, the screening to second order in the impurity potential is not negligible at all, but on the average as large 
as 15 ± 10%. As a guide for the eye we gave the average of Z scr for the three potentials as a dotted line. 

For security we evaluated an alternative expression for Z scr , given by 



kk' 



£fe' [(efc ~ £fe') 2 + a 



212 ' 



(23) 



This expression follows if Eq. H15[l is modified such, that the dipole operator r commutes with the force operator f 1 (t) 
instead of with the statistical operator. While Z scr given by Eq. l|17f) reduces to Sham's expression after taking the 
T — > limit, such a proof is not available for Z scr given by Eq. J22J. On the other hand, the two expressions i|17|) 
and 123|) are equivalent, being related to each other through a partial integration for the derivatives with respect to 
the k and k' vectors. The alternative for Eq. (|18f) becomes 



Z scr = 



37T 2 TO Jq 



dei 



dey 



n k -n k > (£fc-£fc') 2 -a : 



(24) 



If one takes proper care of the higher sensitivity of the expression l|24[) to the choice of the infinitisimal parameter a 
the results turn out to be the same. It can be taken relatively small, much smaller than a typical value of 0.01 for the 
inverse transport relaxation time. In fact, it is the mesh of the integration that determines the lower limit of a. For 
the expressions (|18|) and H24|) it was never larger than 0.00015 and 0.005 respectively. On the other hand, for a = 0.01 
the screening represented by the curves in FIG. ^ reduces by at most 2%. 



IV. REDUCTION OF LINEAR-RESPONSE FORMULA FOR THE SCREENING 



The evaluation of the linear-response formula (jl5|l to all orders in the impurity potential can be achieved by 
restricting the evaluation to a system with one impurity in a jellium, which is in accordance with earlier work by 



others^ In that case f 1 = i[p, h] = i[p, ho + v 1 }. By writing Eq. (|15fl in terms of eigenstates of h 
by q and q' , one can carry out the time integral and finds 



h 



labeled 



Z«^ = l^2 <q \[ T ,n{h)-n(h ) 

qq' 



w> 



t a 



Pq'q 
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= I 53 < - n{h )]\q' >•(-!+ ) = 



3 

qq' 



e<j=<V 



-|ir{[r,n(/i)-7i(/i )]-p} + | 53 < fflMCO ~ n ( h oW > -Pg'g- (25) 

qq' 



The first term in the last line reduces to — Zi because i tr{[r v , n{h) — n(ho)]p^ 1 } 
i tr{{n(h) — n(ho)) [p^, r u }} = tr{n(h) — n(ho)}6 u ^ = ZiS^.^. One arrives at 



i 



Z scr = -Zi + Z cm with Z cm = - J2 < ^1 [r, n(h) - n{h )]\q' > -p q , q . (26) 

qq' 

In view of Eq. @ this would imply a correction term Z COII — Z<j to the cancellation of the bare direct valency Zi. 
The step of subtracting and adding an ia term in the numerator in the second line of Eq. 1251) may look somewhat 
artificial, and the ia factor creates the impression to lead to a zero result in the a — ► limit. But in following Sorbello 
it is seen in the third line that for the e q — e q i terms the ia factor cancels, and the remaining terms give a finite 
contribution^ 

Sorbellc 3 starts from a result obtained by Rimbey and Sorbelloi^ through an evaluation of Eq. J7J and finds after 
some rewritings for Z& 
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Z<j = Imtr 

3nm 



{p 2 (G(6 F )-G°(e F ))}. (27) 



The single particle Green's function G(e) for one impurity in a jcllium, with h = ho + v 1 , and the free electron Green's 
function G°(e) are given by 

GU) = and G°(e) = — — with h = (28) 

w e + ia - h v ' e + ia- h 2m v ' 

This form for Z^ is Sorbello's Eq. (12) and it implies a cancellation of Zi present in his Eq. (7). In order to distinguish 
our result for Z& from Sorbello's Z<\ we keep the notation Z coiI . 



V. COMPARISON OF THE TWO DESCRIPTIONS 

While the description by Rimbey and Sorbello is rather involved and the result l|26[) is obtained in a few lines, it 
is worth while to compare the final expressions. We first evaluate Z covr to lowest order in the impurity potential, for 
which we take a screened Coulomb potential, Eq. (|21l) . By using Eq. (|16[) and the equality 



r, n(h) 



-L I ds n{h)e sh pe- sh (1 - 71(h)) , (29) 
in which the electron mass m = ^ in atomic units, one finds straightforwardly for Z® OII 

ZLr = 5 53 < n(h) - n(h )]\k > k = -An^v kk = Z, = Zj(pot), (30) 



A- 



in which a quantity Zi(pot) = — Ait^/epv kk is defined to be used below. For the screened Coulomb potential this 
quantity is equal to Zi = 1, but this is not the case for other potentials. This result from an explicit calculation 
follows also if one writes the sum over the free space states \k > as a trace and uses the equality given in the sentence 
just below Eq. (J2SJ. 

Similarly one finds for Sorbello's Z^ to lowest order in the impurity potential, writing the trace in Eq. I|27|l in terms 
of free space states labeled by k, 

Z T = -^-Im / d 3 k k 2 G° k (e F )v kk G° k (e F ) = -iiry/^v kk = Z t . (31) 
onrri J 

This is obtained by using the following two equalities, 

A G o (e) d 1 = G o {e)G o {e) and l] m l m G° k (e F ) = -7r6(e F -e k ). (32) 

oe k oe k e — e k + ia a->o 
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Apparently, to lowest order in the impurity potential the two final expressions Z CQTr and are equal and they 
reproduce the bare valency of the migrating ion. Although the complete expressions are not equal, an 'almost' equality 
can be derived. We rewrite Z C orr by applying Eq. Ij29(l both for h and ho. After inserting a complete set of free electron 
states in the /iq term in Eq. (|26|) and carrying out the integral over s, one finds 



7 £ w 

^corr — 

3 m 

qq' 



— 53 - n q )Pqq' ~ 53 nk " ^ ~ nfc ") < > k " < fc "l <? ' > ) ' p 9'9 

gg' fe" 

— 53 (lmG 5 (e F )p w - - 53lmG^(e F ) < #" > k" < fe"]?' > ) ■ p q > q 

qq' k" 
e q =e<t' 

— Im 53 < g|p (G(e F ) - G°(e F )) |g' > -p q , q . (33) 



99' 

In the transition from the first to the second line the T — > limit was taken, for which (3n q (l — n q ) = —-§^-n q — > 

5(e q — e F ) = — — ImG g (eF)- Both Sham and Sorbello give their elaborated expressions in this T — > limit. The 
similarity of this last line with Sorbello's Zj, Eq. (J37J, is striking. The factor of 2 reflects whether the electron spin 
degeneracy has been accounted for explicitly or not. In fact, if in the last line of Eq. I|33l) the states q and q' are 
replaced by the unperturbed ones k and k' , it reduces to Sorbello's expression. This implies an intriguing equality 
indeed, and it shows that the two descriptions are closely related. 

VI. THE CORRECTION TERM IN TERMS OF PHASE SHIFTS 

For the evaluation of the correction term Z co „ as it is defined in Eq. (|26|l one needs the states \q>. These states 
are the eigenstates of a system with one impurity in free space. It is known that the scattering states \ipk > for this 
system, which have a one-to-one correspondence to the free space states \k >, are exact solutions of the Schrodinger 
for one impurity in free space as well. It appears that the evaluation of Z cmT becomes relatively simple if one uses 
the scattering states instead of the true eigenstates. We return to this point below. 

The expansion of the scattering state < r|^fe >= i/'fe( r ) m terms of spherical harmonics is given by 



4-7T 

^« = ^^ e Yl(k)Ri(r,k)Y L (f). (34) 



The angular momentum label L combines the labels t and m, so L = £m, and Rg(r,k) is the radial solution of 
the Schrodinger equation at the energy for a spherically symmetric potential v centered at the origin. For r 
outside the range of the potential Rg(r,k) can be written in terms of the scattering t matrix tt — — r sin<^ exp(i6t) 
as je(kr) — iktehf (kr) , where St are the phase shifts. This means that for a plane wave Re(r, k) — > je(kr). The box 
normalization in the system volume 17 induces a discrete set of k values. In the properties to be presented below a 
delta function normalization will be used, which means that in Eq. I|34|l the system volume fl has to be replaced by 
8ir 3 . Using the expansion Eq. the equality < k'\v\ipk >= tk'k which holds for scattering states \q >^ \ipj~ >, 

and the overlap property for scattering states 

<k'\A >= 6(k ~ k>) V Y L Ck')Y*(k) (l - iktt) , (35) 



one finds for Z c , 



I ^dk(-S(e k -e F ))F( 



2 / , , 2 3 d 



Z ™t = -^zz: I k 3 dk(—S(ek-e F ))F{ek) = z [F{e F ) + -e*—F{e F )), (36) 



in which the function F(ek) is given by 

F(e k ) = \ 53(^ + 1) ( sin 2S e + sm2S e+1 ^j ^ cos 2 {S e - 8 t +i) + lj , (37) 

and F(ek) = hF( e k)- Crucial steps of the derivation of Eq. I|36|l are given in Appendix lAl This expression can 
be evaluated simply, because it is just a function of the phase shifts of the impurity potential at the Fermi energy. 
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•Zcorr as it is given by Eq. (|36[1 has to be compared with the lowest order expression, obtained by the replacements 
| q >— ► \k > and \q' >— ► \k' >. For the sake of a proper comparison this expression has to be evaluated in a similar 
way, by the use of scattering states. This way one obtains 

z° coir = k3dk {i~ k 6iek - e ^) F °^ = ~ + r^^ 0(eF) ) ' (38) 

in which the function F°(ek) is given by 

F°(e k ) = ~^(2£+l)sm2<^ = ^(2£+ l)sin^cos^, (39) 

i i 

and F°(e k ) = lF°(e k ). The right-hand sides of Eqs. (EHJ and ^ can have different numerical values, because 
the two elaborations are different in character. Comparison of these numerical values can be considered as a test of 
the error made in using scattering states instead of the true eigenstates. Another test of this error is obtained by 
evaluating the simplified lowest order result for Z corT explicitly, to be denoted as Z°® TI , using the expansion jTSfr. One 
finds 

Z°°„ = tr{n(h) - n(h )} = - YYW + 1) smS e cosfc = -F°(e k ) ^ - V(2^ + l)S t = Z F , (40) 

7T * — ' 7T 7T ' — ' 

I i 

in which Z F stands for the Friedel sum. The inequality in Eq. (|40|l must be attributed to the use of scattering states 
instead of the true eigenstates. The difference between these two types of states has been stressed by FentonSi and 
commented on by the present autho»2£. A scattering state is prepared such that it is an unperturbed state at t = — oo, 
which develops under the influence of a scattering potential located at some position in the system. An eigenstate 
has to be constructed using the boundary conditions of the system in addition to the properties of the potential. In 
terms of scattering theory eigenstates contain the influence of backscattering by the boundaries in addition to the 
information about the scattering by the potential. It is clear that the third member of Eq. I|4U|I reduces to Zp in the 
small phase-shifts limit Se — > 0. 

Sorbello's equation for Z^, Eq. (|27|l . can be evaluated through the use of scattering states as well. We first give 
the form corresponding to the eigenstates of h, denoted as usual by \q >, followed by the result obtained by using 
scattering states. 

z d = -^imVV<#'>t' 2 G!..M<t'H9>G,M 

q ft' 

-ht.C e|<fc * <(ep - e * )#0(Et) - lii-A^ <41) 

As above, for scattering states the potential matrix element becomes equal to the t matrix element t k r k . The energies 
e k and e q are equal, being connected by a delta function as shown in Eq. (|35f) . Interestingly, this rewritten Zd is 
equal to our Z® OII given by Eq. (|38|l . This again shows the close relationship between the results obtained through 
the simplified approach presented here and Sorbello's results. 

Sorbello calculated Zd of Eq. (|27|l using the square well potentials described in TOII To that end he derives his 
rewritten form Eq. (17), which we reproduce in a slightly different notation as 

Z d = ^ f° r 2 dry j (2i + l){Rj(r,k ¥ )~jj{k F r)) 
3tt Jo ^ 

4ep \— * , „ n 9Se , Ak-pVnrn , 

i 

In the derivation has been used, that ^ = h-v = h-e F + e F -v, G-G° = GvG°, {h-e F )G = -1, — ImTr{G— G } = 
£ £ (2£+ l)ffJ £F , ImG°(r,r) = -fc F El JlW, and ImG( r, r) = : - k F J2 L R 2 L (r), with R L (r) = R £ (r,k F )Y L (f). For 
r < tq the radial solution Ri(r, k F ) oc ji(k v r), with k v = yfef + «o ■ For reasons of a proper comparison we evaluated 
Zd according to Eq. I|42|l up to t — 2, because Sorbello restricted himself to £ max = 0. 
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VII. NUMERICAL RESULTS 



The expressions obtained will be evaluated for the same square well model potentials as were used in £11111 in that 
employing a slight generalization of the potentials used by Sorbello 3 . Results for Zd/Zp and Z COII /Z® otI are shown in 
FIG. |21 Because the Friedel sum for the model square well potentials is rarely equal to unity, the use of ratios gives 
the proper measure for the screening, in which we follow Sorbello. 

It is seen that for the stronger potentials, with ro = 2/A, Sorbello's boldly dashed curve lies somewhat lower than 
the dashed curve for the present description, while they display almost equal results for the higher kp values. The 
curves for the weaker potentials lie lower than those for ro = 2/A. Although they differ considerably for smaller kp 
values, the curves approach each other for higher kp. Both results imply a decrease of the inaccuracy related to the 
\q >^ \tpk > replacement in the present description for states with increasing k values. This is reasonable, because 
larger k values correspond to smaller wave lengths, which probe the scattering potential more precisely, while the 
boundary effect decreases. As a guide for the eye the average for the two well widths are drawn as dotted lines. 
From the present description one comes to a direct valency of 0.85 ± 0.15 on the average, while this is 0.91 ± 0.10 for 
Sorbello's description. So it appears that the \q >^ \ijjk > replacement is not too crude in determining a measure for 
the amount of screening. The screening mentioned by Sorbello is based mainly on the ro = 1/A potentials, because 
for these potentials his restriction to € max = is reasonable. If we correct for the higher I values we find 0.82 for 
Zd/Z-p in aluminum instead of his 0.75. He used the latter value in mentioning a screening of 25%. For the sake 
of completeness we remark, that if one would compare the boldly dashed curve with the values given in Sorbello's 
Table II one would observe considerable differences. This is due to the fact that for the ro = 2/A potentials the t = 1 
and 2 terms contribute significantly. Taking everything together the available models and descriptions end up at a 
screening between 5 and 30%. Comparing with the largely metallic-density dependent result of Ishida, covering the 
entire range of no screening to complete screening, the present result can be considered as rather conclusive, in that 
complete screening is excluded^ 

In order to get some more insight in the kp dependences shown in FIG. we display the quantities Z c olT , Zi(pot) 
and Zp in FIG. 01 It is seen that all curves have a positive slope, apart from the one for Zp for the weaker ro = 1/A 
potential. This is certainly related to the seemingly deviant Zd/Zp curve in FIG.|3|and the fact that the Zd/Zp 
curves have been obtained by an accurate numerical evaluation of Eq. (|42|l . while the Z COTT /Z® m . r results 'suffer' from 
the \q >— > \tpk > replacement. The effect of this replacement is larger for the weaker potential, which can be seen 
from both Z COIT /Z® OII curves in FIG. Eland the Z® OII curves in FIG. 21 This is understandable if one realizes, that the 
effect of the potential on the wave functions increases with its strength, while the effect of the boundaries remains 
unchanged. Further we observe that Z® ml , Zi(pot) and Zp lead to different curves. The difference between Z® OI1 . and 
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FIG. 4: Zcovr, Zi(pot) and Zf plotted as a function of Icf, for the two square well potentials. 

Z,(pot) as shown in FIG.^is certainly due to the \q >— > \tpk > replacement. However, Eq. (|3t)|l implies that for the 
screened Coulomb potential one should find Z,(pot) = Zj? orr = Zj = Zp = 1. This means that a difference between 
Zi(pot) and Zf uncovers some limitation of the use of the model square well potential. Although the differences in 
shape of the three potentials are known, in FIG. [S] we show their shapes for a certain fcp value, for which we chose a 
value in the middle of 0.7. 
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VIII. THE FORCE EXPRESSION OF BOSVIEUX AND FRIEDEL 



The starting expression of Bosvieux and Friedel for the driving force is 

F T=0 = - < *|V Rl (7 ei + SV)\* >, 



(43) 



in which the state |\E' > is a solution of the Schrodinger equation for the system in the presence of an applied field. 
This means that one has to solve the time dependent Schrodinger equation 



=nimt) = ( H + S v(t))*(t). 



(44) 



The subscript T = is added by the present author in order to distinguish this force from the force given in Eq. J2J • 
The derivation of Eq. (0), in which a solution of the Liouville equation is used, has been given in the literature many 
times. Because the approach through the system wave function is typical for the theory of Bosvieux and Friedel, and 
because some questions can be raised about their solution, we give crucial steps of the derivation. We solve Eq. ll4^|) 
by using the interaction representation for ^(t), defined as 



*j(t) 



AHt 



*(*)• 



The equation for ^j(i) becomes 



(45) 



(46) 



After integrating this equation and using that for t — ► — oo the system is in the ground state of the unperturbed 
system Hamiltionian H, one finds for ^(t) linearly in 5V 



*(*) 



-iHt 



dt'e lHt ' : 



<5^(t')e- lfft '*/(-oo) + e^ t * / (-oo) 



(47) 



With SV(t) = 8Ve at , applying the substitution t — t' = s, and considering an arbitrary time in the present, so t = 0, 
this becomes 



/■OO 

*(0) = * = -i / dte-^ H+a ^5Ve iH Hi(-oo) + *j(-oo). 



(48) 



If one calculates matrix elements with this >, the factor e _l£;o00 in the state oo) > drops out so that just 

the ground state jV'o > oi H remains. 
By this one finds for Eq. (|4l5|l 



poo 

F T =o =-< ^o|V Rl ((5F)|Vo >+i dt e~ at < ^ \(V Rl V ei )e- iHt SVe im \i; > + c.c. 

Jo 



ZidE + i / dt e- at < Vol (V Rl V e , 



iHt 



SVe 



■Hi 



Ho > 



(49) 



Because H commutes with the coordinates Rq, in SV only the electron coordinates survive. Using the definition of 
the force operator, Eq. (|SJ), and the hermitian property of H it is clear that the time dependence can be applied to 
F op as well, and Eq. I|49(l can be written as 



F T=0 = Z t eE -ie dt e~ at < ifo] [F op (i), E • ^ |Vo > 



(50) 



Interestingly, Eq. I|5U|I is precisely the zero-temperature equivalent of Eq. J2| . This becomes even more clear if one 
writes down the form which shows up after the reduction of Eq. I|5U|) to single particle states denoted by \q >. 



F T =o = ZteE - ie 



dte~ at J2 <q\U 1 (t),E-r 



Q> 



(51) 
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The force operator f 1 is defined in Eq. J5J. At T — the sum over the single particle states has a sharp cut-off at 
e q = ep ■ The finite temperature equivalent of Eq. H51|l can be written as 



Z,eE 



le 



dt e- at tr\n(h) f^E-r 



(52) 



in which the Fermi-Dirac distribution n(e) has been inserted, see Eq. (|13|) . Clearly, Eq. (|52l) is completely equivalent 
to Eq. (|12f> of the present text. By this electromigration theory can be considered as being unified. Apparently, the 
starting formula of Bosvieux and Friedel was correct, but these authors did not recognize its precise contents. In fact, 
they wrote down surface-integral terms, by this not appreciating the hcrmitian property of the system Hamiltonian. 
This property implies that these terms are zero, but their full-screening results were derived from these terms. In 
addition, they missed the power of their starting formula, Eq. H43J) ■ by taking the a — > limit in too early a stage of 
the derivation. 



IX. CONCLUDING REMARKS AND PERSPECTIVES 



The amount of screening of the direct force on a proton in an electric-current carrying metal has been shown to lie 
between 5 and 25%. By this the full-screening prediction of Bosvieux and Friedel has been invalidated, completely in 
agreement with an earlier result obtained by Sorbelloi^ On top of that, the surface integral terms used by Bosvieux 
and Friedel to derive their full-screening result appear to be zero, due to the hermitian property of the Hamiltonian. 
Interestingly, it has been shown explicitly that the starting expression of Bosvieux and Friedel for the driving force is 
the zero-temperature limit of all linear-response expressions used in the literature since their introduction by Kumar 
and Sorbello£ 

All calculations up to now use a jellium model for a metal, or are not applicable to transition metals^ In view of 
the description presented it becomes feasible to account for real metallic effects. These effects have been accounted for 
for the wind force to a large detail, but for the direct force this was too much involved up to now£ Such a development 
would be interesting, because this may lead to an explanation of a measured result which has not been understood 
yet, up to now. For most hydrides a direct valency for the hydrogen has been measured which is of the order of unity. 
However, in Nb(H) a direct valency was found of about 0.44. Such a deviating value may arise from multiple scattering 
effects of the electrons around a proton surrounded by metallic atoms, which can be accounted for in a finite-cluster 
description. It is worth while to investigate this possibility, because in the development of the description of the 
wind force surprising positive values for the wind valence in V(H) and Nb(H) were found, which were in agreement 
with experiments^. The surprise comes from the fact, that in a system composed of a finite cluster embedded in 
a jellium the electron dispersion relation is still free-electron like, from which one would expect a negative wind 
valence. The calculated result must be due to the rather strong multiple scattering effects, which were accounted for 
explicitly. Similarly, it is worth while to develop a finite cluster description for the direct valency, which is a straight 
generalisation of the impurity in a jellium description implemented so far. This is a feasible development if one uses 
the simplified treatment presented above, of which it has been shown that the expressions can be evaluated in terms 
of the scattering phase shifts of the constituent atoms. A slightly more explicit hint can be found elsewhere^ 

Acknowledgement I want to thank Jacques Friedel for an interesting conversation and an extended correspondence 
over the recent months. 



APPENDIX A: CRUCIAL STEPS IN THE DERIVATION OF EQ. 

For the evaluation of Z COIT defined in Eq. (|26l) one needs the momentum matrix element p q ' q and the matrix element 
of the commutator with the statistical distributions. For p q i q one writes 

P«'o = IT? E^'^'^W / d 3 rRl,(r)pR L (r), (Al) 

in which Eq. I|34|l with £1 — > 8ir 3 has been used for the wave functions and Rl{y) = Rt(?, fe)Y^(r). If one represents 
the scattering potential by a square well with depth vo the inner radial solution is a Bessel function as well, so 
Ri(r, k) = Aiji(k v r), with k v = \/k 2 + vq, and one finds 

J Y£,(f)p R L {r)df = i t '- i k v n L , L R t u{r, k), (A2) 



13 



in which the equality 



Y L *,(f)p j L (v)df = i" ~"k T> L , L j t {kr) 



has been used and 



D 



L'L 



dkY£,{k)kY L {k). 



(A3) 



(A4) 



The double I- label in Eq. (|A2(I refers to the fact, that the factor Ag is not changed by the momentum operation and 
the angular integration, so that 



A 



A e , 

If one substitutes Eq. I|A2(1 in Eq. (|A1|) one obtains 



k v Rf((r, k) = k v —R e >(r, k) — > k {ji>{kr) - ikt e hj,(kr)) 



P 9 ' 9 = -Y, Y L'(k r )Yl(k)-D L , L I^f(k',k) 

7T * " 



LL' 



with 



7^f(fc',fc)= / r 2 drk r R},(r 7 k r )R e , e (r,k), 



and in which k r — k v inside the range of the potential and k r — k — kp outside of it. Using the equality 

7T(5(fc - k') 



r 2 drji(k'r)ji(kr) = 



2k 2 



it appears to be possible to reduce the integral I RR (k' , k) to 

T RR (h' h\ - klT S ( k ~ T RR 
1 tU \ K i K > — p i 

in which 

Ipf 1 = 1 — iktg + ikt\, — 2iktfikt\, . 
Now we turn to the other matrix element in Eq. (|26|) . Using Eq. (|16fl one finds 



< q\ [r, n{h) - n{h )]\q' >= - ds < q\ 



r,n(h)e sh ve' sho (l~n{h )) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



We remind the reader that in this equation h and v just refer to the system with one impurity, so to h 1 and v . Now 
we use the equality l|29(l and the following related equality 



i 

m 



ds p(s)e 



Oh 



(A12) 



twice, one time for h and one time for ho. By that the commutator in the right hand side of Eq. IjAllj) can be written 
as 



r,n(/i)e s,l ue- ;5 ' l0 (l - n(h )) 



= -L n (h)( f ' ds'p(s')(l - n(h)) 
ds'p(s')e sh ve~ sho + e sh ve- sho p (s - pn(h ))^j (l - n(/i )) ■ 



e ah ve -sho 



(A13) 



It will be clear that p(s) in the first and second term refers to h. 

Now we develop the qq' matrix element of this operator as it occurs in Eq. (| All|> . In that we will make use of the 
property proven above through Eq. I|A6|) with l|A9|) . namely that the energies e q and e q r in the matrix element p qq r 
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are equal, and of the equality of the energies e k and e q in the overlap < k\q >, see Eq. (135 j l . In the first and second 
term we have to insert two complete sets, one q set \q"' >< q"'\ and one k set \k" >< k"\. In the third and fourth 
term one needs the complete k set only. This way one writes for the qq' matrix element in the left hand side of Eq. 
(lATTl) 



i0 

< q\[r,n(h) - n(h )]\q >= n q 



Xm p 99'"(^ _? V")ife"fc'" + /Z t *k"k k "(l -n k ") 



{l-n k „)<k"\q' >, (A14) 

1 J ' L 

qtll k n 



in which the potential matrix element v qk i — v k , ql between an exact scattered state \q >= \i)) k >, see Eq. (|34H . and 
an unperturbed state, a plane wave \k' >, has been replaced by the corresponding t matrix element t k , k . 

By now we have developed the means for bringing Z COTI in a manageable form. One has to take the inner product 
of the matrix element given by Eq. i|A14fl with ^p q >q and to carry out the summations. Because the summations are 
equivalent to integrals and the absolute values of all fc-vectors involved arc equal through delta functions, one just 
has to carry out the angular integrations. We write for Z COII in Eq. iffi^! 

Zcorr = \J & Jdk' Jdk" (J dk"^j < q^nih 1 ) - n(h )]\q' > -p q , q , (A15) 

in which the right hand side of Eq. I|A14|) is supposed to have been substituted. The angular integration over k'" 
applies to the first term in the right hand side of Eq. (|A14J| only. The product of statistical factors which shows up 
can be simplified and in the T — > limit be written as follows 



I3 2 n k (l - n*)(^ - n k ) = lpJL(n k (\ _ „ fc )) _> \-^5{e k - e F ). (A16) 



If one substitutes this equality, uses Eq. l|A6Jl with Eqs. (|A9j) and (|A10(I . and Eq. 1A14|I . accounts for the factor of 2 
due to the spin degeneracy and carries out all angular integrations, one finds for Z COII of Eq. (|A15(I 

Z m rr = -z / k 2 dk—— 5(e k - e F ) V" D LL / • T> L > L (1 - ikt e >) x 

3irm J Q de k ^ 

{t%I%* + t$) eiff = --i- ^ k 3 dk(^-S(e k - e F ))F(e k ), (A17) 
37rm J \ de k > 

in which has been used that 

I** = cos(5/ - 6t>) e - l{5 '- 5 >'\ (A18) 

and that 

E D "ll>D£, l - ^ Ul + 1) S e4+1 + i h>,l-A ■ (A19) 

mm' 

In the right hand side of Eq. (|A17J) one recognizes the second member of Eq. Ij36j) ■ by which the derivation has 
been completed. The derivations leading to Eqs. 138|) . I|4l)[) and (|41|l are similar, and they are simpler as well. 
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